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Abstract 



The PHENIX experiment uses semileptonic decay channels to measure open and closed heavy flavor cross sections across the 
rapidity range -2.2 < y < 2.4. High luminosity data are now available for p+p, d+Au, Cu+Cu and Au+Au collisions at -\fs~^=200 
GeV, and for Au+Au collisions at -\/s NN = 62 GeV. We discuss recent d+Au results for open heavy flavor and J/if/ production, and 
discuss their implications for the cold nuclear matter contributions to heavy flavor production in heavy ion collisions. 
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1. Introduction 

At low Bjorken momentum fraction (x) in nuclei, where gluon densities are very high, saturation effects are 
expected to become important UJ. The modification of gluon densities in high energy collisions involving nuclear 
targets is inherently interesting because it contains information about such gluon saturation effects, and also because 
it determines the initial conditions in high energy nucleus-nucleus collisions, where the hot matter effects can be 
understood only after the initial conditions are known. 

Charm and bottom quarks are attractive as probes of the hot matter created in nucleus-nucleus collisions because 
their large mass prevents them from being created thermally in the hot medium. They are created only in hard 
processes that occur during the nuclear crossing, which is short in high energy collisions due to the large Lorentz 
contraction of the colliding nuclei in the collision frame. At RHIC, for example, the crossing time is ~ 0.13 fm/c, 
which is smaller than the thermalization time of the QGP 

The use of heavy quarks and quarkonia to probe the hot matter created in nuclear collisions would be straight- 
forward if heavy flavor production scaled with the number of nucleon-nucleon collisions, because then the initial 
distributions could be obtained from measurements in p + p collisions. However there are a number of effects in p+A 
collisions that modify the initial distributions of open heavy flavor and quarkonia, including shadowing. These are 
called cold nuclear matter (CNM) effects. They have to be studied using p+A data (or at RHIC, for technical reasons, 
d+A data). For this reason, PHENIX has carried out a program of measurements that includes so far p + p, d+Au, 
Cu+Cu, Au+Au and, very recently, Cu+Au collisions. 

2. d+Au results 

PHENIX has recently released new results on open heavy flavor (HF) production in d+Au collisions measured via 
semileptonic decays at mid rapidity |2|. The pj dependence of the nuclear modification is shown in Fig.[T]for central 
and peripheral collisions. Enhancement is observed for the 0-20% most central collisions in the pj range 1-5 GeV/c, 
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while little or no enhancement is observed for peripheral collisions. At high pj, where strong suppression is observed 
for HF electrons from Au+Au collisions, the measured RdAu ls consistent with one. Thus the suppression observed in 
Au+Au collisions can be attributed, within the uncertainties on RdAu, to hot nuclear matter effects. 
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Figure 1. The RdAu for HF electrons measured at -^s^n = 200 GeV for a) the most central collisions and for b) the most peripheral collisions 0. 

Published PHENIX d+Au J/if/ data from the 2008 RHIC run span the rapidity range -2.2 to 2.4 in 12 rapidity 
bins, and have been divided into 4 centrality bins Q. Fig. [2] shows the J /if/ RdAu as a function of rapidity for each 
of the four centrality bins. The forward rapidity data for the most central collisions show suppression of about 
50%. The centrality dependence of the forward rapidity data reveals an interesting feature, which is illustrated in 
Fig. [3] [3 1. In this plot the centrality integrated RdAu on the horizontal axis represents the overall modification. The 
ratio of the modification for the most central and most peripheral collisions, Rep, on the vertical axis characterizes 
the centrality slope of the modification. The points show the locations of the PHENIX measurements at twelve 
rapidities, with the systematic uncertainties represented by ellipses. The lines on the plot are loci generated from a 
Glauber model of the d+Au. collisions in which the modification for each nucleon-nucleus collision is assumed to 
have certain simple mathematical dependencies on the nuclear thickness Ta(tt) at the nucleon impact parameter rj. 
The three mathematical forms shown are exponential, linear and quadratic. The plot shows that the forward rapidity 
modifications require a dependence on Ta(jt) that is quadratic, or greater. 

PHENIX has recently released measurements of the pj dependence of the J/tfr RdAu as a function of centrality 2). 
The data are compared in [4] with various model calculations. A comparison of the centrality integrated data with two 
models is shown here in Fig. [4] One calculation [5|, shown as the dot-dashed line, uses a model of the color dipole 
breakup based on a parameterization of HERA data, a Cronin broadening that is parameterized from low energy data, 
and a shadowing correction obtained from the nDSG nuclear parton distribution function (nPDF) parameterization. 
The calculation over-predicts the suppression, but otherwise reproduces the shape of the pj distribution at forward 
and mid rapidity. However it does a poor job of reproducing the shape of the distribution at backward rapidity. In the 
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Fi gure 2. y/i/' RdAu versus rapidity for four centrality bins l3l . 



45 second calculation (6] |3 shadowing is parameterized using various nPDF sets, and the effect of collisions of the cc 

46 dipole with nucleons during the nuclear crossing is represented by an effective breakup cross section. The example 

47 shown in Fig.|4]as the dashed line uses the nDSG nPDF set, with an effective breakup cross section of 4.2 mb. There 

48 is no added Cronin enhancement. This second calculation describes the data reasonably well at low pi at forward and 

49 mid rapidity, although it under-predicts the data at higher pj. However at backward rapidity this calculation also fails, 
so At backward rapidity and low pj (where x ~ 0.1 for the parton in the Au nucleus) production occurs in the 

51 anti-shadowing region, while at high p T (x ~ 0.3) production moves towards the EMC region. No direct evidence 

52 of an EMC effect has been reported in the gluon distributions, and few constraints on the gluon distributions exist in 

53 this region. Furthermore there is large disagreement in the modification of the gluon density between nPDFs. The 

54 nDSg nPDF set includes no suppression in the EMC region, and only a small anti-shadowing effect, while the EKS98 

55 nPDF exhibits a suppression similar to that observed in the quark distributions, and a larger anti-shadowing effect. 

56 The lack of a strong anti-shadowing effect combined with the absence of an EMC effect in the nDSg nPDF causes the 

57 calculated /? ( ; A „ to remain roughly constant with increasing p T . For the EKS98 nPDF set, the larger anti-shadowing 

58 combined with the inclusion of an EMC effect causes a decrease in the calculated RdAu as pj (and correspondingly 

59 x) increases p4J, worsening agreement with the data. The explanation of the pj shape at backward rapidity remains 
eo unclear. The inclusion of a Cronin effect in [5 1 that seems adequate to describe the more forward rapidity data does 
61 not allow a good description of the backward rapidity data, suggesting that the disagreement may be due to the nPDF 
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Figure 3. RdAu plotted versus Rep (!3j. The data points show the location for the 12 rapidity bins measured by PHENIX, the ellipses represent the 
systematic uncertainties. See the text for discussion. 

sets at high x, or possibly to some other physics effect not considered. 

The J /if/ pj broadening relative to p + p collisions can be characterized by plotting the difference of the values of 
(Pj) between the d+Au and p + p distributions. This difference is shown as a function of centrality, characterized by 
(N co n), in Fig. [5] The dependence on centrality is very similar at forward and backward rapidity, with a hint that it is 
larger at mid rapidity. However it should be noted that (pSj) is about 20% larger at mid rapidity than it is at forward or 
backward rapidity for both p + p and d+Au collisions [4j, so that the difference plotted in Fig.[5]roughly scales with 
the size of {p\). 

3. Heavy Ion Results 

There are new semileptonic decay open heavy flavor PHENIX results for heavy ions that include Cu+Cu Raa data 
at both mid and forward rapidity. 

In Fig. [6] and Fig. [7] the mid rapidity Cu+Cu Raa data are compared as a function of N co n with data from d+Au 
and Au+Au collisions. The data are shown for two p T ranges. Where they overlap, the three data sets are found to 
display similar behavior with N co u within uncertainties. It is noteworthy that the the d+Au data and the peripheral 
Cu+Cu data have very similar modifications, including some enhancement at N co u ~ 10. 

At forward rapidity (1 .4 < |y| < 1 .9) , recent results on the Raa for semileptonic decays of heavy quarks to muons 
in Cu+Cu collisions have now been released JH). The Raa for the most central collisions is shown in Fig. [8] The 
suppression in central collisions is found to be stronger than that seen at comparable N part at mid rapidity for Au+Au 
collisions, suggesting the possibility that CNM effects are larger at forward rapidity. The observed suppression is 
consistent with a calculation |9| that includes the effects of heavy-quark energy loss and in-medium heavy meson 
dissociation, as well as cold nuclear matter effects due to shadowing and initial state energy loss due to multiple 
scattering of incoming partons before they interact to form the cc pair. 

PHENIX has recently released measurements of the Jj^i Raa at V^aw = 62 and 39 GeV in the rapidity range 
1.2 < \y\ < 2.2 IfTOl . These are compared in Fig.|9]with the existing 200 GeV data [11|. The modifications are similar 
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Figure 4. The J/ij/ Rjau versus pj at backward, mid and forward rapidity |4|. The theory curves are described in the text. 



at the three energies, but it should be noted that CNM effects are expected to be quite different for J /if/ production at 
the three energies lfl2l . We do not yet have comparison d+Au data at 62 and 39 GeV. 

Taken together, the ALICE Pb+Pb data at ^fs^ = 2.76 TeV (BJ, the NA50 Pb+Pb data at ^fs^ = 17.3 
GeV [14], and the PHENIX Au+Au data at 39, 62 and 200 GeV provide a J/i// data set that spans a large energy 
density range and a broad range of initial conditions and charm production cross sections. At sufficiently high cc 
pair production rates, a significant fraction of the J /if/ yield is expected to be due to coalescence of c and c quarks 
from different hard processes in the same nuclear collision, and comparison of data at LHC and RHIC energies will 
be necessary to constrain this contribution. Once the necessary CNM reference data are obtained for the PHENIX 39 
and 62 GeV data, and the ALICE 2.76 TeV data, the available data should provide strong constraints on theories of 
J / if/ modification in heavy ion collisions. 



4. Summary 



Recently released PHENIX RdAu data at 200 GeV for HF electrons at mid rapidity show enhancement for central 
collisions in the range pj = 1-5 GeV/c, with no evidence of significant modification at higher pj. These results 
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Figure 5. The enhancement of (p|) for d+Au collisions over p + p collisions, plotted versus collision centrality at backward, mid and forward 
rapidity (4). 



ga indicate that the strong HF electron suppression observed for Au+Au collisions at medium to high pj is due to hot 

99 matter effects. 

100 New results have now been released for the centrality and pj dependence of RdAu f° r J/ ft production at 200 

101 GeV. The data were measured in rapidity intervals -2.2 < y < — 1.2, |y| < 0.35 and 1.2 < y < 2.2. Together with 

102 previously published RdAu vs centrality data at twelve rapidities, these complete the PHENIX J/t[r R^au results from 

103 the 2008 RHIC run. Two notable observations have been made about the data. First, the centrality dependence of the 

104 modification at the most forward rapidities is found to be stronger than linear in the nuclear thickness. Second, the 

105 pi dependence of RdAu at backward rapidity is not described by available calculations, suggesting perhaps a problem 
ioo with the nPDF parameterizations for x values in the anti-shadowing and EMC region. 

107 There are new PHENIX preliminary HF electron Raa data from 200 GeV Cu+Cu collisions. Where they overlap, 

108 d+Au, Cu+Cu and Au+Au show very similar behavior with N a ,n, including enhancement of the d+Au and Cu+Cu 

109 data at N co u ~ 10, at lower pj. Newly released HF muon Raa data from 200 GeV Cu+Cu collisions show suppression 
no that is stronger than is seen at mid rapidity at the same N part , perhaps consistent with stronger CNM effects at forward 
in rapidity. 

112 PHENIX has now released J/ift Raa data from 39 and 62 GeV Au+Au collisions. The modifications are found to 

113 be similar to that at 200 GeV. However CNM effects are known to be strongly energy dependent, and CNM reference 
iu data will be needed before conclusions can be drawn about the differences between hot matter effects at these energies. 
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Figure 6. Heavy flavor electron modifications at mid rapidity for 
Cu+Cu collisions compared with those for d+Au and Au+Au. In 
all cases the modifications are integrated over the range pr = 1 - 3 
GeV/c. 



Figure 7. Heavy flavor electron modifications at mid rapidity for 
Cu+Cu collisions compared with those for d+Au and Au+Au. In 
all cases the modifications are integrated over the range pj = 3 - 5 
GeV/c. 
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Figure 8. The pj dependence of the Raa for semileptonic decays of heavy quarks to muons in Cu+Cu collisions at 1.4 < |y| < 1.9 Q5). The theory 
curve |9| is discussed in the text. 
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Figure 9. Comparison of newly released Au+Au J/ifi Raa data at 62 and 39 GeV with existing data at 200 GeV. 



134 
135 



GeV, Phys.Rev. C86 (2012) 024909. 
[9] R. Sharma, I. Vitev, B.-W. Zhang, Light-cone wave function approach to open heavy flavor dynamics in QCD matter, Phys.Rev. C80 (2009) 

054902. 

[10] A. Adare, et al., J /if/ suppression at forward rapidity in Au+Au collisions at ijs NN =39 and 62.4 Ge\ arXiv: 1208 . 2251 

[1 1] A. Adare, et al., J /{ft suppression at forward rapidity in Au+Au collisions at ysjv/v = 200 GeV, Phys.Rev. C84 (201 1) 054912. 

[12] C. Lourenco, R. Vogt, H. K. Wohri, Energy dependence of J /if/ absorption in proton-nucleus collisions, JHEP 02 (2009) 014. 

[13] B. Abelev, et al., J/psi production at low transverse momentum in Pb-Pb collisions at sqrt(sNN) = 2.76 TeV, Phys.Rev.Lett. 109 (2012) 

072301. 

[14] R. Arnaldi, J/psi production in p-A and A-A collisions at fixed target experiments, Nucl.Phys. A830 (2009) 345C-352C. 



8 



